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Abstract 

Although diagrammatic perturbation theory fails for the dynamical-mean field theory of the 
double-exchange model, the theory is nevertheless ^-derivable and hence thermodynamically con- 
sistent, meaning that the same thermodynamic properties are obtained from either the partition 
function or the Green's function. We verify this consistency by evaluating the magnetic suscepti- 
bility and Curie temperature for any Hund's coupling. 
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The dynamical mean-field theory (DMFT) formulated in the late 1980's by Miiller- 
Hartmann [1] and Metzner and Vollhardt Q] has developed into one of the most power- 
ful many-body techniques for studying electronic models such as the Hubbard 0, Q| and 
double-exchange (DE) models. This theory is believed to become exact in the 

limit of infinite dimensions and to capture the physics of correlated electron systems even in 
three dimensions. Recent work on dilute magnetic semiconductors has used DMFT to study 
variants of the DE model 



lOl with less than one local moment per site. In this paper, 
we reach the surprising conclusion that, unlike for the DMFT of the Hubbard model ^, a 
diagrammatic perturbation theory containing only electronic degrees of freedom fails for the 
DMFT of the DE model. Nevertheless, we show that the theory remains ^-derivable in a 
more restrictive sense, which still implies that the partition function and Green's function 
produce consistent results for thermodynamic properties such as the magnetic susceptibility 
and Curie temperature. 

The Hamiltonian of the DE model is given by 

H = -t^2 (c\ a c ja + c] a c ia ) - 2 J H 51 s, ■ Si (1) 

<j,i) * 

where c\ a and Q a are the creation and destruction operators for an electron with spin a 
at site i, s, = (l/2)c| a cr Q/ gc i/3 is the electronic spin, and Sj = Sirii is the spin of the local 
moment. Repeated spin indices are summed. Within DMFT, the effective action on site 
above T c in zero field is given by 

A eff (m) = -T^2c 0a (w n ){G (w n )~ 1 5 al 3 + Jn^ap • m}c 0f5 (iis n ), (2) 

n 

where Jh = JnS, v n = (2n+ 1)ttT, co a {iv n ) and co a (iu n ) are now anticommuting Grassman 
variables, and G^iiv^) is the bare Green's function containing dynamical information about 
the hopping of electrons from other sites onto site 0. 

Because A e s(m) is quadratic in the Grassman variables, the full Green's function G(iu n ) a p 
at site may be readily solved by integrating over the Grassman variables, with the para- 
magnetic result 

G(i Vn )L = ({G„(i,„)-'/ + J s a ■ m }~ l ) m = fry /, (3) 

where J is the unity matrix in 2 x 2 spin space. The average over the orientations m 
of the local moment is generally given by (C(m)) m = / dQ m P(m)C(m), where P(m) oc 
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FIG. 1: (a) The bare vertex function; (b) and (c) Compact diagrams that contribute to <I> for the 
electronic effective action A' efi on the right with their associated self-energies on the left. 

Tr^exp(— y4 c fj(m)^ is the probability for the local moment to point in the m direction. 
Above Tc, P(m) = 1/Att is constant. Consequently, the paramagnetic self-energy is given 



by E(iz/ n ) = G Q (iu n ) 1 



G{iv ri 



JftGo^iVn). Expanded in powers of Jh and G{iu n ), we 



find 



1 
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+ Jl = J^G(iu n ) - j^G(iv n f + 2J^G(tv n ) 5 + 

DO 



2G(iv n ) V 4 G(^n) 2 
On a Bethe lattice, these relations are closed by the analytic expression 

W 2 



G (iv n ) = z n I - —G(iu n ), 
lb 



(4) 



(5) 



where z n = iv n + fi and W is the full bandwidth of the non-interacting, semicircular density- 
of-states. We denote the full spin dependence for later use. 

Diagrammatic perturbation theory is customarily formulated in terms of the bare vertex 
function (l,n;m) l3a ' 5fi sketched in Fig. 1(a) with tu m = 2rmrT. The bare vertex function 



may be associated with the two-particle interaction in the purely electronic effective action 



may 





A' eS = -T^2coa{iu n )G (iu n ) x c® a {iv n 



rp3 

— y E co<*(iv n + iu m )c p(iv n )r {0 \l,n]m)P a ' M c 0K {m)c $(iv l + iu m ). (6) 

l,n,m 

Hence, the bare vertex function must satisfy the crossing symmetries T^(l,l + m;n — 
l) Sa 'P K = r^(n +m,n;l - n)^ K;fo = -r(°'(l,n;m) WB . There are two ways to calcu- 
late r(°)(z, n; m)' ' . First, we can take the Jh — > limit of the full irreducible vertex 
07,n: m) l3a ' &K obtained from the Bethe-Salpeter equation for the magnetic susceptibility 
[a, Alternatively, we can associate the lowest-order, J H contribution to the parti- 

tion function Z = (Tr(exp(— A eff (m)))) m with the contribution to the partition function 
Z' = Tr(exp(— A' eS )^j, sketched as the compact diagram in Fig. 1(b) (with internal lines given 
by the bare Green's functions Go(iv n ) a /3). Both methods yield the same result: 

r(°)(/, n- mf a ' 6K = \[3Jl{<Jpa ■ (T SK 8 m ,o ~ vsa ■ <rp K 6 ln }, (7) 

which satisfies the crossing symmetries. 

However, replacing A e ff(m) by A' eS produces an inequivalent theory Q]. For example, 
expanding Z and Z' in powers of Jh yields the results 

Z = Z {l - J H E G 0(^n) 2 + ^J H E G 0(^) 2 G (^n) 2 + O(JS)}, (8) 

Z' = Z {l - J H E G {iv n ) 2 + |j H E GviivtfGSvnf + 0(J£j}, (9) 

which disagree to order J H . Hence, it is not possible by averaging over the local moments 
to reduce the Hund's coupling to an effective two-particle interaction. In other words, 
the Hund's coupling produces fourth and higher-order electronic interactions that require 
higher-order vertex functions in the electronic action. 

A theory is usually said to be ^-derivable if a functional $({G(«i/ n )}), constructed from 
the sum of compact diagrams in terms of the full Green's functions and the bare vertex 
functions, can be found to satisfy the condition S(z^„) aj g = 5$/<5G(iz/ n ) aj g. As discussed by 
Baym [15|, a ^-derivable theory may readily be shown to be thermodynamically consistent, 
meaning that thermodynamic properties can be evaluated either from the Green's function 



or from the partition function Z. For a ^-derivable theory, the partition function Z or free 
energy — T log Z may be constructed in terms of $ from the relation 

- logZ = $ - ]TTr{s(*z/ n )G(^n)} + ^Trlog{G(^„)}, (10) 

n n 

which is stationary under variations of G_{iv n ). Whereas Baym's original work was intended 
for systems of interacting Fermions and Bosons, the notion of $-derivability has been ex- 
tended to systems of interacting electrons and spins [l(| and to disordered alloys [ftI ]. 

From the discussion above, it is clear that even if it exists, $ cannot be constructed 
in terms of the bare vertex functions. When the action contains only two-particle inter- 
actions such as for the Hubbard model, then the first two terms in $ are represented by 
the compact diagrams on the right-hand side of Figs. 1(b) and (c) with the correspond- 
ing self-energies T;(ii/ n ) aa = 5$ / '5G(iv n ) aa sketched on the left-hand side. Not surpris- 
ingly, substituting our earlier expression for the bare vertex function produces the cor- 
rect first-order self-energy E W(iz/ n ) = J^G(iiy n ) but the wrong second-order self-energy 
E'®{iis n ) = -{J^/3){2G{iiy n )J2 l G(iiy l ) 2 + G(iu n f). Notice from Eq.® that the correct 
second-order self-energy S^(w„) = —J^G(iu n ) 3 does not involve a Matsubara summation. 
Hence, the DMFT of the DE model is not ^-derivable in the strict diagrammatic sense 
stated above. 

Despite the failure of a diagrammatic expansion in powers of E (0) , a func- 
tional $({G(ii/ n )}) can still be constructed to satisfy the condition Y,(ii> n ) aa = 

/ 5G{iv n ) aa . Starting from Eq.(j3J) and Dyson's equation for the self-energy, we find that 
8Yi(ivi) aa / 5G{iv n ) = (K~ 1 )fP + 5i n 5 a pG(ii> n )~ 2 , where K_ is the Jacobian 

with a n = G (w n )~ 2 — Jn and b n = Go(* z/ n) -2 — <^h/ 3 - This Jacobian can be inverted with 
the general result 

dG(iv n )pp 3b n [ 2a n - 3b n J 3 - 2 £ r l/6 r kK v y 

It can be shown [13] that the right-hand side equals —TT(l,n;m = 

QjaatfP where 

r(Z, n; m) l3a ' ,SK is the full irreducible vertex of the Bethe-Salpeter equation. The func- 
tional $ must exist because the curl of the self-energy vanishes: 5T l (iu n ) aa / 8G(iui) pp — 
5E(iui)pp/SG(iu n ) aa = 0. 



By construction, $W (second order in Jh) is represented by the compact diagram in 
Fig. 1(b) and is given in terms of the bare vertex function by 

= -|E r( ° ) (^;0)^G(^)aaG'(^)/ 3/3 = ~ J -f {E G (^)aa 

x (G(iu n ) aa - Gr(ii/„)aa) - E G ( iu n)l a ~ 2 E G{iv n ) aa G(iv n )^ }, (13) 

n n ^ 

where a is the opposite spin to a. After expanding and integrating Eq. (fT2|) [l8^ . we find 
that $*- 2 ) (fourth order in Jh) is given by 

y ^ 4 n n Z,n,r 

2 1 

x (G{iu n ) aa -G{w n )ao) + rE G K)L( G ( w i)«« - G(ivi)aa) \. (14) 

Unlike $^ 2 - ) cannot be represented by a compact diagram involving only the bare vertex 
functions. So far, all of our results are valid for any lattice topology. 

We have verified the thermodynamic consistency of the DMFT by calculating the mag- 
netic susceptibility from both the Green's function and the partition function. With a 
magnetic field H = Hz coupled to both the local moments and the electrons, the effective 
action becomes 

A cS (m) = -T^c 0a (iv n ){G (iv n )~p + ( J H m + -Hz) ■ (T a p\cQp(iv n ) - (3HSm z . (15) 

Parameterizing the bare inverse Green's function as G (iz/ n ) _1 = (z n + R n )L + QnQ. z an d 
using Eq.© for the full Green's function, we solve for R n and Q n on a Bethe lattice from 
the expression 

Rnl + Q n a z = -^r({(*n + R n )L + (j H m + (Q n + H/2)z) ■ a.y 1 ) m (16) 
To linear order in the field, R n and Q n satisfy the implicit relations 

R n = Zn + Rn _„ (17) 

16 (z n + R n y- Jg' 1 ' 

n H(z n + R n ) — 2J R Mi m R n H 

Wn " 2(z n + 2R n )U n 2' [ ^ 

where 

32J| Rl 

3WHz n + R n )(z n + 2R n y 1 ] 



After integrating exp(— A e g(m)) over the Grassman variables, we find that the probability 
for the local moment to point along m is 



P(m) oc expj^logj 1 - ^—^^t^j^ ) + /3HSm z \ oc exp((3J eff M lm m z 

{ n \ \Z n + R n ) ~JuJ ) 



(20) 



where the local-moment order parameter M; m = (m z ) is solved from the condition M; m = 
J e ffMi m f3/3. The electronic order parameter M e i = 2(sq z ) is obtained from the summation 
M el = -(32T/W 2 )J2 n Qn- The total susceptibility is then given by the zero-field limit of 
X = (SMi m + M e i/2)/H. To calculate the susceptibility from the partition function, we 
first expand Z to second order in H and M/ m and then use x — (T / H)d\ogZ / dH\u=o- 
The latter technique is formally equivalent to evaluating the susceptibility from the Bethe- 
Salpeter equation 

These two sets of calculations do indeed produce the same magnetic susceptibility, which 
may be written as 

*= h r- (jffew + ^ (Gl(r) " Ga(r)) + ^ Gl(r) - (21) 

3 f T 

SAT) =S+ 2#r(GiCO - G 3 (T)), (22) 
where the functions G\(T) and G^iT') are formally given by the Matsubara sums 

32 i? 2 
Ci(T) = -y E (2n + + 2R n )U n ' (23) 

G 2 (T) = — — ^2 __}^r>__ ^24) 

The Curie temperature T c is solved from the condition Gi(T c ) = (W/Ju) 2 . 

Previous results in the Jh — > oo limit are reproduced [3| by taking /i = sgn(p — 
1) J H + Sfi where < W/(2y/2) and p is the electron filling (p — 1 means one electron 
per site). The general expression for the magnetic susceptibility shall be studied in a future 
publication. We pause here to note that the effective spin S e s(T) may be either larger 
or smaller than S depending on the sign of the Hund's coupling Jh- The temperature 
dependence of S e g(T) and the deviation of 1 — ( Jn/W) 2 Gi(T) from T — Tq are both caused 
by electronic correlations that are absent in a local- moment system y|. The second and 
third sets of terms in Eq. (|21|) correspond to the Pauli susceptibility of the electrons. 

Although $({G(w n )}) has no simple diagrammatic interpretation, the existence of this 
functional means that we may still use Eq. (jlO|) to establish the thermodynamic consistency 



of the DMFT of the DE model. Diagrammatics may be recovered for a more sophisticated 
model where the classical local moments are replaced by fully quantum-mechanical operators 
and we introduce an additional propagator corresponding to those local spins. It may also 
be possible to develop a more complex diagrammatics for classical local spins in terms of 
higher-order vertex functions. 

Finally, we note that whereas any conserving theory (in the sense of Baym and Kadanoff 
is thermodynamically consistent, it is not true that all thermodynamically consistent 
theories are conserving. Indeed, that is the case here since the DMFT violates the Ward 
identities associated with charge and spin conservation. 
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